Draft version July 9, 20 1 2 

Preprint typeset using I4TgX style emulateapj v. 5/2/1 1 



CLUSTERING PROPERTIES OF BzK- S ELECTED GALAXIES IN GOODS-N: ENVIRONMENTAL QUENCHING AND 

TRIGGERING OF STAR FORMATION AT Z - 2 

Lihwai Lin Mark Dickinson 2 , Hung-Yu Jian \ A. I. Merson 4 , C. M. Baugh 4 , Douglas Scott 5 , Sebastien Foucaud 6 , 
Wei-Hao Wang ', Chi-Hung Yan Hao-Jing Yan 7 , Yi-Wen Cheng 8 , Yicheng Guo 9 , John Helly 4 , Franz Kirsten 10 , 
David C. Koo ", Claudia del P. Lagos 4 , Nicole Meger 5 , Hugo Messias 12 , Alexandra Pope 9 , Luc Simard 1314 , 

Norman A. Grogin 15 , Shiang-Yu Wang 1 

Draft version July 9, 2012 

ABSTRACT 

Using a sample of Z?z^-selected galaxies at z ~ 2 identified from the CFHTAVIRCAM near-infrared survey of 
GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), 
and stellar mass (M*), and the clustering of galaxies as a function of these parameters. For star-forming 
galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR oc M* a with 
a = 0.74 ± 0.20 down to M* ~ 10 9 M Q , indicating a weak dependence on the stellar mass of the star formation 
rate efficiency, namely, SSFR. We also measure the angular correlation function and hence infer the correlation 
length for sBzK galaxies as a function of A/*, SFR, and SSFR, as well as K-bwA apparent magnitude. We show 
that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring 
the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies 
ranges from 4 to 20 h~ l Mpc, being a strong function of M%, M*, and SFR. On the other hand, the clustering 
dependence on SSFR changes abruptly at 2 x 10~ 9 yr" 1 , which is the typical value for "main sequence" star- 
forming galaxies at z ~ 2. We show that the correlation length reaches a minimum at this characteristic value, 
and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from 
the predictions of the semi-analytical models. Our results suggest that there are two types of environmental 
effects at work at z ~ 2. Stronger clustering for relatively quiescent galaxies implies that the environment 
has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with 
elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers 
in dense environments. 

Subject headings: galaxiesxlustering - galaxies:evolution - galaxies: high-redshift - large-scale structure of 
Universe 



1. INTRODUCTION 

When exactly galaxies formed, and how their evolutionary 
histories are associated with their environment, are among 
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the most important problems in extra-galactic astronomy. 
Recent studies have suggested that the critical epoch for 
building up galaxy mass and shaping galaxy properties is 
around z ~ 2, as the peak i n the cosmic star formation rate 
densi t y lies at this redshift ([Dickinson et al.l l2003t IMo pkins 
120041 Hopkins & Beacoml 120061: ICharv & Popell2010l) . Sim- 
ilar trends have also been discovered for AGN/QSO activity 
dRichards et al.l 12006). Furthermore, the number and stellar 
mass density of quiescent gal axies have been fo und to build 
up rapidly since this epoch (lArnouts et a l. 2007). Therefore, 
it is essential to probe the properties and abundance of galax- 
ies at this epoch and beyond in order to understand galaxy 
evolution. One method which selects z > 1 galaxies is the so- 
called BzK color selection. Specifically at 1.4 < z < 2.5, this 
method has been shown to efficiently separate star-formation 
dominated sourc es from those alrea dy in a passive phase of 
galaxy evolution (Daddi et al. 2004|). 

While large numbers of z ~ 2 galaxies have been routinely 
discovered by deep imaging surveys, we still lack a complete 
picture of the connection between different galaxy properties 
at this era, and in particular the interplay between the host 
dark matter halos and galaxies. One method to quantify the 
masses of dark matter halos that host the galaxies is to mea- 
sure the amplitude of galaxy clustering. Within the cold dark 
matter (CDM) model, it is well known that more massive dark 
matter haloes are more s trongly clustered (e.g., Bau gh etafl 
U999tlMo & Whitd l2002). Despite this, the details of the spa- 
tial distribution of galaxies may be complex due to the addi- 
tional physical processes affecting galaxies during their for- 
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mation and evolutionary histories, for instance, super-novae 
and AGN feedback (e g .,| kav et al]|2002HBenson et aH l200l 
ISomerviile et al. 2008|). However, under simple assumptions 
or through halo occupation distribution (HOD) modeling, 
one can infer the h ost halo mass of galaxies through their 
clustering strength dBertin et alJ 120021 : iRravtsov et al.l I2004t 
iZheng et alJl2005l:lTinker et al.ll2005h . 

Observations at z < 1.5 suggested that the clustering 
amplitude depends strongly on galax y properties such as 
morphology, colo r , and luminosity (|Norberg et al.l 120021 : 
iZehavi et alJlIOOH iCoil et all 12001 Ide la Torre et alJl201lT) . 
At z ~ 3 and higher redshifts it has also been found that 
the clustering strength increases with UV-continuum lu- 
minosity and hence SFR for UV-selected Lyman Break 
Galaxies dGiavalisco & Dickinsonl200lUFoucaud et all2003l: 
l Adelberger et alJ )2005t iKashikawa et al.1 120061: iLeeetall 
l2006t lOuchi et alj|2004t lYoshida et al.ll2008l) . However, the 
situation for z ~ 2 populations is less clear. There have been 
several attempts to study the clustering of Bz/f-selected pop- 
ulations. Most of these focused on the dependence of clus- 
tering on /f-band magnitu de and galaxy type s (star-forming 
BzKs vs passive BzKs: iKong et al.l 120061 lHavashi et all 
200^lBlanc et alJl2008l:lHartlev et al.ll2008l:lMcCracken et alJ 
2010). Only a few studies have probed the clustering depen- 
dence on other properties such as stellar mass (M„), star for- 
mation rate (SFR), and spe cific star formation rate (SSFR = 
SFR/ s tellar mass) at z ~ 2 dFoucaud et al.ll2010t [Wake et alJ 
1201 UlMagliocchetti et aT]|201 UlSavov et al.11201 II) . This was 
partly due to the difficulty in obtaining deep and yet wide-field 
near-infrared (NIR) imaging, which is critical for the BzK se- 
lection. The samples were typically drawn from either small 
deep fields or shallow wide surveys, and therefore lacked the 
dynamic range required for sampling the properties of inter- 
est. 

Thanks to the new-generation of wide-field NIR imagers, 
deep NIR galaxy surveys have now become possible, allow- 
ing one to probe both bright and faint populations of galaxies 
simultaneously. In this work, we carry out the first system- 
atic study of clustering as a function of various galaxy prop- 
erties for Bz/f-selected galaxies down to K s ~ 24.0 (AB mag), 
selected from the multiwavelength data i n the Great Ob- 
serva tories Origins Deep Survey (GOODS: Gia valisco et alJ 
120041) North field (GOODS-N). The relation between clus- 
tering strength and SSFR is particularly interesting because 
we can gain insight into the relation between the halo mass, 
and thereby environment, and star formation efficiency. Pre- 
vious work suggests that environment plays an important role 
at l ow redsh i fts in quenching the star formation of galax- 
ies dDressleri 119801: iBalogh et alJ fl998t iGomez et all 120031: 
Kauffman n et alj|2004f). but that this e ffect reverses at z ~ 1 
(Elba z et all 120071: ICooper et alJ [2008h . However, a recent 
study bv lOuadri et al.l d2012l) investigating the quiescent frac- 
tion as a function of local density found in contrast that the 
SFR-density relation continues at least out to z ~ 1.8. The 
role of environment at z > lis thus still under debate, although 
different results may be attributed to different definitions of 
the 'SFR-density' relation (see §4.2). Instead of using lo- 
cal density measurements, in this work we use the cluster- 
ing strength as a probe of halo mass and thus environment. 
This has the advantage of avoiding noisy density measure- 
ments due to either photometric redshift errors, or incomplete 
sampling in the spectroscopic redshift samples given that the 
method relies on the projected distributions alone. 

The paper is organized as follow. The data sets, BzK sam- 



ples, and the methods used in this paper to compute the stellar 
mass and star formation rate are described in §2. The results 
of the relation between SFR, SSFR, and M*, and the cluster- 
ing analysis are presented in §3. We discuss the implications 
of our results in §4, followed by a brief summary in §5. 

Throughout this paper we adopt the following cosmology: 
H = 100/z km s" 1 Mpc" 1 , Q m = 0.3 and fl A = 0.7. We adopt 
the Hubble constant h = 0.7 when calculating rest-frame mag- 
nitudes. We use a Salpeter IMF when deriving stellar masses 
and star formation rates. All magnitudes are given in the AB 
system. 

2. DATA, SAMPLE SELECTIONS, AND METHODS 

GOODS-North is one of the most heavily studied extra- 
galactic fields, with rich multiwavelength data sets. The NIR 
observations in GOODS-N were carried out with the Wide- 
Field Near Infrared Camera (WIRCAM) on the CFHT dur- 
ing 2006-2009. These include 27.4 hr of integration in J 
band obtained by a Taiwanese program (PI: L. Lin) and 31.9 
hr of integration in K s band obtained by Hawaiian (PI: L. 
Cowie) and Canadian (PI: L. Simard) programs. Results from 
the WIRCAM A s -band imaging were recently published by 
IWang et alJ d20Tol) ; here we use our own reductions (L. Lin et 
al., in preparation). The data were first pre-proc essed using 
the SIMPLE Imaging and Mosaicking PipeLinE (Wang et al. 
2010), and then combined to pro duce deep s tacks with the 
AstrOmatic softwa re 'SCAMP' dBertinl [2006h and 'SWarp' 
dBertin et al . 2002) PI The resulting 5a limiting magnitudes 
using 2-arcsec diameter circular apertures, reach J = 24.6 and 
K s = 24.0 in the central 420 arcmin 2 . The optical data used in 
this work comes from two sources: one is the HST/ ACS B435, 
V606, '775' an d Z850 v2.0 catalog from the GOODS HST Trea- 
sury Program (Giavalisco et alj|2004l) : the other includes the 
ground-based U band imaging obtained on the KPNO May- 
all 4-m telescope and the BVRiz bands taken with Subaru 
dCapaketalJ 120041) . The GOODS-N has also been imaged 
with IRAC at 3.6, 4.5, 5.8, and 8.0^m , taken as pail of the 
GOODS Spitzer Legacy program (PI: M. Dickinson). The 
area used in this work is limited to the ACS coverage of 10' 
x 16' = 160 arcmin 2 . 

The photometry for each band is done using a soft- 
ware package wit h object template-fitting method (TFIT; 
Laidler et al. 2007) based on the ACS z-band detections. For 
each object, TFIT constructs a template using the spatial po- 
sition and morphology of the object in the ACS z-band image. 
Such a template is then convolved with the PSF of other low- 
resolution images and then fit to the images of the object in 
other bands. The best-fitting fluxes are then considered as the 
final fluxes of the object in low-resolution bands. The ACS 
photometry is measured in dual-mode with SExtractor while 
photometry in other bands is measured by TFIT with z- band 
template (N. Grogin, et al. in preparation.) 

Our star-forming BzK galaxies (hereafter sBzKs) are se- 
lected using HST/ACS B and z bands, and the CFHT 
WIRCAM K ? band , following the BzK method developed by 
iDaddi et all ( 120041) . We add 0.3 mag to the z-K color of the 
dividing lines that separate the sBzK galaxies, passive BzK 
galaxies (hereafter pBzKs), stars, and foreground/background 
galaxies (hereafter non-Bz/f s), in order to account for the pho- 
tometry offsets due to different filter systems and zeropoint 
calibrations. The sBzK galaxies are selected to be those ob- 

16 http://www.astromatic.net/ 
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jects with 



(z-K)-(B-z)>0.1, 



(1) 



while pBzK galaxies are defined as those objects with 

(z-K)-(B-z)< 0.1 n(z-K)> 2.8. (2) 



Stars are also identified by the relation 

(z-K)-(B-z) < -0.2+(B-z) x 0.3. 



(3) 



The B-z versus z-K distribution of our sample is shown 
in FigureQ] We require signal-to-noise ratios S/N > 5 in both 
the z and K bands Moreover, we further apply a magni- 
tude cut of [4.5] < 25.3 mag, corresponding to a 5-cr limiting 
mag in [4.5] so we can obtain reliable stellar mass measure- 
ment. For objects undetected in the B-band, if their BiaUm-z 
and z-K colors satisfy the pBzK criteria, we classify them as 
pBzK galaxies, where B2 a u m is the B band 2a limiting mag- 
nitude. On the other hand, we treat the remaining objects 
without detections in B as unclassified BzK galaxies (here- 
after uBzKs). In total, we classified 4496 sBzKs, 45 pBzKs, 
179 uBzKs, 341 stars, and 5791 non-BzKs. The BzK sample 
used in this work probes galaxies fainter than previous clus- 
tering studies of BzKs by 0.5-1 mag in the /T-band. 

The photometric redshifts (hereafter z p hot) were obtained by 
fitting ground-based UBVRizJK photo metry using th e pub- 
lic "BPZ" photometric redshift code (lBenitezl l2000). The 
rms of our z p hoi when compared to the spectroscopic redshifts 
of the same objects is about 0.07(1 +z) with 6.5% outliers 
at 0.1 < z < 3.0. Spectroscopic redshifts used for our z p hot 
calibration were compiled fr om several p u blishe d and unpub- 
lished s ources, parti c ularly iCohen et al.l (120001). Iw"irth et al.l 
d2004l) . iReddv et al.l d2006bl) . iBarger et all (120081) . and D. 
Stern et al. (in preparation). This spectroscopic redshift cat- 
alog consists of 4308 objects up to z ~ 6.7, among which 
639 objects satisfy our sBzK selection. Although the spec- 
troscopic sample contains objects with A'-band mag down to 
25.0, they are systematically brighter than the full BzK sam- 
ple used in this work. The K-b&nd magnitude of the 639 
sBz/f-selected galaxies with spectroscopic redshift measure- 
ment peaks around 23.0, while that of the full sBzK sample 
peaks at K ~ 24.0. 

The Zphoi distributions of the sBzK and the non-BzK popula- 
tions, together with the spectroscopic redshift distribution of 
sfiz/f-selected galaxies which have a spectroscopic identifica- 
tion, are shown in Figure [2] It can be seen that although the 
original design of the sBzK selection is to pick up galaxies at 
1 .4 < z < 2.5, the redshift distribution of our sample suggests 
that there is non-negligible contaminati on from galaxies a t 
Z < 1 .4 as well as those at z > 2.5 (also see lBarger et al.l20"08h . 
Such contamin ation however does not affect our clustering re- 
sults in Section [3~2l since we take into account redshift distri- 
bution when converting the angular cluster amplitude into the 
real-space clustering strength. For those BzK galaxies which 
do not have a photometric redshift measurement because of 
detection in limited bandpasses or whose redshifts fall out of 
the redshift range 1 < z p hoi < 4, we assign their reshifts to be 
the median Zpi,,,, for a given subsample used in the clustering 
analysis. None of our conclusions change if we restrict our 
sample to only those with photometric redshift measurement. 

The ACS B435 band samples rest-frame wavelengths rang- 
ing from 1200 to 1800 A for the BzK galaxies (1 .4 < z < 2.5), 

17 The typical magnitudes with S/N = 5 in our sample are 27.9 in the B 
band, 27.3 in the z. band, and 24.0 in the K s band. 
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FIG. 1 .— B-z vs. z-K diagram for all objects in GOODS-N. The sBzKs, 
pBzKs, foreground/background galaxies (non-Bz.Ks), and stars are classified 
according to the four regions separated by the solid and dashed lines as 
adopted by Daddi et al. (2004) with additional small adjustments to account 
for the filter differences. Stars represent the pBzK galaxies that are detected 
in the B-band while the arrows denote the objects not detected in the B-band, 
but their Bi a n m —z and z-K colors satisfy the pBzK criteria. 
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FIG. 2. — Photometric redshift distributions for sBzK (blue histogram) and 
non-BzK galaxies (green histogram). The blue curve represents the best fit of 
the Gaussian parameterizing the z p hot distribution for sBzK galaxies. Objects 
with photometric redshift value smaller than 0.1 are not shown. For com- 
parison, we also show the spectroscopic redshift distribution for those sBzK 
galaxies that have a spectroscopic redshift identification (red histogram). 
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allowing us to estimate the rest-frame UV luminosity at 1500 
A. According to iDaddi et al. (2004), the reddening for sBzK 
galaxies correlates well with the observed (B - z) color in the 
following way: 

E(B-V) = Q.25(B-z+0.l) AB . (4) 

One potential caveat of this approach is that the above em- 
pirical relation is derived based on a brighter sBzK sample 
with K < 22 (Dadd i et all 120041) . and has not been directly 
tested for fainter /f-selected galaxies. However, the color- 
reddening relation for sBzK galaxies was define d assuming 
the dust attenuation law from Cal zetti et alj ( 2000). and mid- 
and far-infrared observations (Reddy et al. 2006a, 2012J) have 
shown that this law works well, on average, for typical L* 
Lyman break galaxies at z ~ 2, which significantly overlap 
with fainter s feg-selected samples . Considering the Calzetti 
extinction law (ICalzetti et al.l l2000). we compute the SFR by 
converting the extinct ion corrected 1^(15 00 A) using the fol- 
lowing equation from lDaddi et al.l ( 120041) : 

SFR(M©yr _1 ) = L„(1500)/(8.85 x lO^Hz" 1 ). (5) 

To estimate the stellar mass , we fo llow a similar method to 
that adopted in Da ddi et ail d2004l) . but using the observed 
IRAC 4.5 fim magnitude and [3.6] - [4.5] color instead. We 
parameterize the stellar mass in the form of : 

/^(M,/M )=ll.O-O.4([4.5]-fl 1 (z)) + a 2 (z)([3.6]-[4.5]). 

(6) 

The a i and «2 values given in Table Q] are estimated based on 
a calibration of Eq. [6] against the stellar mass d erived from 
the SED fitting code 'FAST' dKriek et al.1 12009) for a sub- 
set of bright sBzK galaxies, which have decent detections in 
multiple bands and spectroscopic redshift measurements. Al- 
though the depth of the 4.5 /im data allows us to probe the 
BzK sample down to M» ~ 10 9 M©, our sample is also limited 
by the depth in B— , z~, and K— bands which are used for the 
BzK selection. Similar to Eq. [6] one can also estimate stellar 
mass using the K— band mag and the z-K color as adopted by 
Daddi etaLl (120041) . As a result, we are likely missing galax- 
ies with red z-K colors in the stellar mass range 10 Mq< 
M„< 10 10 M Q given our depth in K. In this work, we carry 
out the analysis down to M*~ 10 9 M©, and we will discuss 
the effect caused by the incompleteness. 

3. RESULTS 

3.1. SFR - M* and SSFR - M* relations 

We begin by showing SFR versus stellar mass for a sub- 
set of the BzK sample with 1.8 < z p hot < 2.2 in Figure 
[3] (left panel). Although the SFR is estimated using the 
UV light corrected for dust extinction, rather than using the 
full SED fitting, there is a clear sequence of sBzK galax- 
ies in the SFR-M, relati on, as seen in other work that in- 
clude s IR SFR tracers (IWuyts et al.l 120111: IRodighier o et al.1 
120111) . Using SFR oc M°, the slope a for our sBzK galax- 
ies is found to be 0.74 ± 0.20 over the stellar mass range 
10 9 M©< M»< 10 11 M©. Both the slope and the normal- 
ization are in good agreement with previous studies at z ~ 2 
(IDaddi et al.1 120071: IPannella et al.1 120091: iWuvts et al.l l20TTt 
Rodig hiero et al.ll201 ll) . However, this work extends the stel- 
lar mass limit further down to 10 9 M©. As mentioned in §2, 



it is possible that we are missing galaxies with red z-K col- 
ors, namely, the relatively quiescent galaxies, in the low mass 
end. If this is true, we would overestimate the SFR for less 
massive sBzK galaxies, and therefore the slope of our derived 
SFR-M* relation should be regarded as a lower limit. 

Also shown in Figure [3] is the l-a distribu tion of the 
SFR-M * relation at z ~ 2 as predicted by the [Lagos et ail 
d201 ll) GALFORM model, implement ed on the halo merge r 
trees of the Millennium Simulation (Springel et al. 2005). 
The GALFORM model is a semi-analytica l model of galaxy 
formation that was originally introduced in lCole et al.l (2000. 
see Baugh 2006 for a review). The SFR and M* are directly 
taken from the GALFORM output. The IM F adopted in the 
model is a Kennicutt IMF dKennicuttlll983l) . so we add 0.3 
dex to both the SFR and M* to convert them into a Salpeter 
IMF. It appears that the model predicts larger scatters in SFR 
and under-predicts the SFR for a given stellar mass compared 
to our observations by roughly a factor of 2-4. 

In the right panel of Figure[3]we plot the SSFR as a function 
of the stellar mass. Our data suggest that the sBzK galaxies are 
forming stars with slightly weaker depend ence on stellar mass 
than what is observed at lower redshifts dNoeske et al.ll2007l: 
iCowie & Bareeril2008l) . The slope of our SSFR-M* relation 
(-0.26 ±0.20) is fla tter than other work based on optical to 
mid-IR SFR tracers dFeulner et alj|2005t iReddv et~aTl2 006a: 
lErb et alj|2006t iPapovich et al.ll2006t IRodighiero et al.ll201oi) 
at z ~ 2, but is consistent with the results derived from radio 
SFR tracers dDunne et all 120091: IPannella et al.1 12009). This 
flattening of the SSFR-M* relation at z ~ 2 has recently bee n 
interpreted as the dawn of downsizing (Pan nella et al.1 2009). 
in the sense that star formation had not ceased in massive 
galaxies at these redshifts. 

It is worth noting that some of the sBzK galaxies, in partic- 
ular for more massive ones, have very low SSFR consistent 
with that of pBzK galaxies, as revealed in Figure[3] This may 
be partly due to the contamination of pas sive galaxies that ar e 
being scattered into the sBzK selection (Ou adri et al.| [2007). 
leading to a potential underestimation of the median SSFR 
for intrinsically massive star-forming galaxies. In addition, 
the inclusion of passive galaxies in the sBzK sample may also 
affect the results of cluster ing o f galaxies with low SSFR for 
which we will discuss in 33.21 Nevertheless we expect that 
the contamination rate should be low given the much lower 
number density of passive galaxies compared to that of sBzK 
galaxies and hence our main results regarding the SSFR (or 
SFR) - M* relation and the clustering analysis should not be 
significantly affected. 

3.2. Clustering Properties 

We measure the angular correlation function (ACF) for the 
sBzK galaxies using data and random ca talogs with the esti- 
mator proposed by Landy & Szalayj d 19931) : 

uj(6) = [DD(6)(n R /n D ) 2 -2DR(6)(n R /n D ) + RR(6)]/RR(6), 

(7) 

where DD(6>), DR(6>), and RR(6>) are the number of data-data, 
data-random and random-random pair counts with separations 
between and + 80, and no and n R are number of galaxies 
in the data and random catalogs. Figure [4] shows the ACF 
measured for the GOODS-N sBzK sample, which is binned 
according to /f-band magnitudes, M*, SFR, and SSFR. The 
error bars in each bin are based on the la Poisson statistics 
of the pair counts in the bins and hence are not correlated 
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FIG. 3. — Left: SFR and stellar mass correlation for Bz.K galaxies in GOODS-N. Right: SSFR and stellar mass correlation. sBz.K galaxies are shown using small 
black dots. The green triangles and blue filled circles represent the median values of the star-forming galaxies in the GOODS-N sample and in the GALFORM 
model predictions respectively. Star-forming galaxies in the GALFORM model refer to those galaxies with SSFR greater than lCr" yr - '. pBzK galaxies in 
GOODS-N are shown using red stars. The ones with arrows represent those pBzK galaxies that are undetected in B. The upper limits of SFR or SSFR are 
estimated using the 2-cr B-band limiting magnitudes. The error bars represent the l-cr distribution in each stellar mass bin. 



across the bins. We fit our ACF with a power law 

lj(6)=A uj (0 1 -~«-C), (8) 

where 7 is fixed to be 1.8 and C is the integral constraint which 
accounts fo r the finite region o f the sky probed in the sample. 
Following Roc he et alj ([1999), we estimate the integral con- 
straint using 



C- 



T,RR(9)& 
ZRR(9) 



-0.8 



(9) 



which gives a value of C = 7.86 in our case. The ltr error 
on the amplitude A u is computed from the covariance matrix 
output by the IDL fitting program "mpfitfun". 

To convert the 2-D clustering amplitude A u into the 3- 
D correlation length rp, we use Limb er's inversion dLimberl 
[19531 EEgliocchetti & Maddox 1999), assuming a Gaussian 
redshift distribution for our BzK sample. Ideally, we should 
fit a Gaussian to each subsample to derive the peak and the 
width of its photo-z distribution. However, the sizes of some 
subsamples are too small to derive meaningful fitting results. 
Instead, we set the mean of the Gaussian to be the median 



of the photometric redshift distribution for each subset, and 
adopt a global width of 0.47 derived from the deconvolution 
of the photo-z width of the whole sBzK sample with the pho- 
tometric redshift error (see Fig. |2). However, it is important 
to note that ro is sensitive to the adopted width of redshift dis- 
tribution, being larger for broader width for a given clustering 
amplitude A u . Therefore if the width varies among different 
subsamples, it might introduce bias in the trend we see below. 
To take into account this effect, we also compute ro in the 
case of a width of 0.3 which is a good estimate of the mini- 
mum value of the width across our subsamplj^ and find that 
ro decreases by 22% compared to the case of width = 0.47. 
We then add this 22% error in quadrature when estimating the 
error of ro. 

The masses of the halos hosting the galaxies are estimated 
from ro in three different ways. The first and simplest method 
is to assume that each halo hosts only one galaxy, and the 
selected sample corresponds to a narrow range of halo mass. 

18 This also corresponds to the width in the case where approximately 95% 
of the sBz,K galaxies are confined within the redshift range of 1.4 < z < 2.5. 
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FIG. 4. — The angular correlation function selected in (a) S'-band magnitude, (b) stellar mass, (c) star formation rate, and (d) specific star formation rate for 
sBzK galaxies. In each panel, the colors represent the measurement of samples selected in different bins. The lines represent the best fit power law using a fixed 
slope of 7 = 1.8, with the integral constraint incorporated. Bins with 8 < 0.001 deg are excluded from the fit in order to reduce the contribution from very close 
pairs. For clarity, only four selected bins are plotted in each panel. The fitting results for the full subsamples are give in Table [2] The values of SSFR which 
appear in the legend are in units of 10~ 9 yr~' . 



The second method also assumes one-to-one correspondence 
between halos and galaxies, but the selected samples are al- 
lowed to be hosted by halos above a certain minimum halo 
mass M m ; n . The third method takes into account the Halo Oc- 
cupation Distribution (HOD), which parameterizes the num- 
ber of galaxies as a function of halo mass, allowing the possi- 
bility that halos can host more than one galaxy. Given our 
small sample size in each bin, we are not able to perform 
a full HOD analysis by fitting the ACF with one-halo and 
two-halo terms. Ins tead of fitting the HOD parameters, we 
follow the recipe of Zheng et al. (2007) and model the num- 
ber of galaxies hosted by a dark matter halo with mass M as 
(N(M)) = 1 +M/M U with M x = 20M min for M > M min , where 
Mmin is the minimum halo mass needed to host one galaxy 
and Mi is the mass when a halo hosts two galaxies. We then 
compute the expected ro for a given M m - m and compare that 
to the ro derived in our sample. For a given M m ,„, one can 
compute the effective ro defined as the sum of the ro of all 
halos with M > M min weighted by N(M), divided by the total 
number of halos of with M > M mi „. We note that the HOD 
adopted here may be oversimplified and not optimized for our 



sample. Nevertheless, it provides a good approximation to 
how the halo mass varies for a given ro due to differences in 
the ro-mass conversions. The derived clustering strength, in- 
ferred halo mass from the above three approaches, as well as 
other sample characteristics, are given in Table [2] 

Figure [5] plots the 3-D correlation length as functions of 
K s , M*, SFR and SSFR for sBzK galaxies. First, we confirm 
previous results that the clustering amplitude increases as K 
brightness increases (Figure and|5^). This correlation be- 
tween the K-bwd brightness and the clustering amplitude is 
often interpreted as clustering strength increasing with galaxy 
stellar mass, as seen at lower redshifts. Since we now have an 
improved stellar mass estimator from IRAC fluxes (see Eq. 

we are able to probe the stellar mass dependence of the 
ACF directly (Figure^ and|3J)). It is evident that there exists 
a strong dependence of the ACF on stellar mass down to the 
stellar mass limit of IO^Mq, which is consistent with the find- 
ings by IWake et al.l (1201 lb . A similar trend is also seen as a 
function of the UV-derived star formation rate (Figure HJ; and 
Figure [5};), as reported by ISavoy et aD (1201 ll) . Moreover in 
Figure [SJ}, we find that the pBzK galaxies are more clustered 
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FIG. 5 . — Correlation length ro for sBzK galaxies as a function of (a) K-band magnitude, (b) stellar mass, (c) star formation rate, and (d) specific star formation 
rate. The black filled circles denote the results of sBzK galaxies in our sample. As a comparison, the results of pBzK galaxies are shown in red triangles. The 
SSFR of pBzK galaxies is arbitrary assigned to be 4x 1CT 1 1 yr~' in the lower-right panel. It can be seen that pBzK ga laxies are more stro ngly clustered than most 
sBzK galaxies. The gray shaded areas represent 1-c predictions from the semi-analytical galaxy formation model of Lagos et al. (201 1). 



than sBzK galaxies with similar stellar masses. 

On the other hand, there is not a monotonic dependence of 
clustering strength on the specific star formation rate. The 
correlation length, ro, decreases with increasing SSFR for 
galaxies with SSFR below 2 x 10~ 9 yr" 1 , but increases with 
SSFR for galaxies exceeding this threshold. 

As a comparison to the HOD predictions, we also show the 
predictions for how the clustering strength depends upon the 
K-bsnA, M*, SFR and SSFR values output by the GALFORM 
model. The GALFORM results, along with 1-cr uncertainties, 
are indicated by the grey shaded region in Figure [5] It can be 
seen that the semi-analytical model predicts similar trends for 
K-bwA magnitude and M» dependence, as seen in the obser- 
vations, but this dependence is typically weaker in the model. 
On the other hand, although the GALFORM model predic- 
tions roughly agree with the ro-SFR relation for galaxies with 
SFR < 50 Moyr" 1 , the models have few galaxies with SFRs 
as large as those we infer for the real sBzK galaxies, limiting 
the range over which we can make the comparison. As a con- 
sequence, it is not straightforward to make comparisons in the 
ro-SSFR relation between model and observations due to the 



narrower dynamical range of SSFR in the model, although 
there is marginal hint of the up-turn signature in the model 
predictions. In order to understand whether this discrepancy 
is due to selection effects, such as color selection, or the physi- 
cal treatment in the model, it is necessary to carry out an anal- 
ysis in the model as it was done in this BzK g alaxy sample. 
This w ill be addressed in a forthcoming paper (Merso n et al.l 
l2012h . 

4. DISCUSSION 

4.1. Comparison of ro between sBzK galaxies and other 
dusty star-forming galaxies at z ~ 2 

In this work, we find that for star-forming BzK galaxies at 
z ~ 2, the clustering amplitude decreases with the apparent 
A"- band magnitude, in good agreement w i th previous stud- 
ies (lHavashi et al.l2007tlHartley et al.l2008HMcCracken et al.1 
2010), but our work extends these studies to less luminous 
populations by almost one magnitude. We also observe a pos- 
itive correlation between stellar mass and clustering strength 
for sBzK galaxies. The inferred spatial correlation length is 
~ 4.5 h~ l Mpc for galaxies with stellar mass around 2 x 10 9 
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FIG. 6. — Correlation length ro of sBzK galaxies binned in SFR from this 
study (SFR increases from left to right among blue symbols; see Table 0, 
compared to the measurements of DOGs (green symbols) and SMGs (red 
symbols) from previous studies. 
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FIG. 7. — The relation between rg and specific star formation rate of sBzK 
galaxies for different stellar mass bins: 9.0 < log(M* /Mq) < 9.5 (red cir- 
cles), 9.5 < log(M*/M@) < 10.0 (green triangles), and 10.0 < /og(M»/M Q ) 
< 10.5 (blue squares). It can be seen that the dichotomy shown in the lower- 
left panel of Figure|5]still holds at fixed stellar masses. 

M© and increases to ~ 13.9 h~ l Mpc for galaxies of 1.5 
xl0 n M Q . 

Moreover, the clustering amplitude also depends on SFR, 
being stronger for galaxies with higher SFR. This trend may 
not be surprising given the good correlation between SFR and 
M*. We find that the least actively star-forming systems, with 



1.0 < SFR/(M©yr _1 ) < 5.0 have r ~ 5.3 h" x Mpc, suggest- 
ing that their typical host halos have masses ~ 2.3 xlO 12 
M Q (based on Mo & White 2002); while those with SFR 
> 100 M Q yr _1 are more strongly clustered with ro ~ 14.1 h~ l 
Mpc and are hosted by dark matter halos with masses above 
2.8-4.5 x 10 13 M Q , depending on the actual HOD models. 

Figure [6] compares our clustering measurement of sBzK 
galaxies to those for other dusty, star-forming systems at sim- 
ilar redshifts. The correlation length we find for sBzK galax- 
ies with SFR > 100 M yr _1 , 14.1 ± 3.5 /T 1 Mpc, is in broad 
agreem ent with that of s ubmillim eter galaxies from recent 
studie s (IWeifi et al.1120091: IWilliams etalll20lU iHickox et alJ 
120121) . which obtained a r of 7.7-13.6 /T 1 Mpc. It is also 
consistent with ro = 12. 97^' ^ n ~ l Mpc fr° m a measurement 
of 24 nm. -selected ( F24 > 0.6 mjy) dust-obscured galaxies 
dBrodwin et alJ 120081) . echoing the finding that highly star- 
forming BzKs and brighter du st-obscured galax i es may actu- 
ally be the same populations (Po pe et al.1 12008: Mege r et alJ 
1201 11) . Moreover, the correlation length of our ULIRG- 
like sBzK galaxies is also comparable to that of far-infrared 
sources detected at 100/im and 160^m with the PACS instru- 
ment on Herschel in the GOODS-South field, which have 
typical correlation lengths of 17.3-19 Mpc, or equivalently, 
12.2 - 13.3 h~ l Mpc, assuming h = 0.7 (iMagliocchetti et al.1 
1201 11) . Despite that the aforementioned populations are se- 
lected using different techniques, all these results appear to be 
converging: at z ~ 2, the most rapidly star-forming galaxies 
are strongly clustered, and highly star-forming BzK's may be 
linked to and/or overlapping with very dusty systems, based 
upon the similarities of their correlation lengths. 

4.2. Environment quenching and triggering of SFR ? 

In the local Universe, it is known that star formation is a 
strong functi on of environment, being less active in dense en- 
vironments (|G6miFitalJ |200l |K^^annital]|200l). It 
is also well established that the star for mation rate is corre- 
lated with stellar mass at all redshifts dFeulner et all l2005t 
iNoeske et all 120071: iPerez-Gonzalez et al.1 120081) . Since the 
stellar mass distributions of galaxies in different environments 
differ, then whether the SFR-environment relation is a purely 
environment effect, or is governed by the stellar mass is still 
under debate. Recent work at lower redshift (z < 1) has sug- 
gested that both mass and environment are responsible for 
shaping th e properties of galaxies and that the i r effects are 
separ able jBamford e t al. 2008; Pe ng et al.ll2010HSobral et all 
1201 11) . The net effect is that in dense environments, or in mas- 
sive dark-matter halos Q3, the averaged star formation rate is 
smaller than that of galaxies located in under-dense regions 
(or hosted by less-massive halos). 

The SFR-environment relation can be probed in vari- 
ous ways, including the SFR-density, SSFR-density, and 
color-density relations. The first two relations can be dervied 
by measuring the averaged SFR and SSFR as a function of 
overdensity, while the third relation usually refers to either 
the relationship between the galaxy color and the overdensity, 
or the fraction of quiescent galaxies as a function of overden- 
sity. As different probes may result in different results, it is 
crucial to specify which method is being quoted wh en inter- 
preting the results as discussed in iPatel et al.l (1201 ll) . For in- 
stance, at z ~ 1, it has been shown that, when considering the 

9 The local density roughly scales with the dark matter halo masses, and 
therefore for the rest of the discussion we will use environment and dark 
matter halos interchangeably 
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populations of both star-forming and quiescent galaxies alto- 
gether, the SSFR-density and color -density relations follow a 
similar trend as that found locally (ICooper et al.ll2007ll2008l) . 
while the SFR-density relation is reverse d in the sense that 
the averaged SFR increases with density dElbaz et alj|2007t 
ICooper eta l. 2008), primarily due to a population of bright, 
blue galaxies, as well as dusty LIRGs in overdense environ- 
ments at this epoc h (ICooper et all 120061: lElbaz et all 120071: 
Coope r"etaill20 08 ) Pi Nevertheless these studies used spec- 
troscopic redshift samples which normally suffer from incom- 
plete sampling of galaxies, leading to large unc ertainties in 
the env ironment measurements. Recent work by Ouad ri et al.1 
d2012l) . who study the fraction of quiescent galaxies using a 
photometric redshift sample drawn from the UKIDSS Ultra- 
Deep Survey (UDS; O. Almaini 201 1, in preparation), claims 
that the color-density relation persists out to z ~ 2 at all stel- 
lar masses, although they caution about the large uncertainties 
due to the errors in the photometric redshift. 

Here we try to address the question of when environment 
comes to play with a different approach from a local den- 
sity measurement: we measure the clustering strength as a 
function of SSFR among star forming galaxies, and also that 
for different galaxy types (star-forming vs passive). If galax- 
ies with lower SSFR tend to located in the denser environ- 
ments, we should see a larger correlation length for galaxies 
with lower SSFR, based on the assumption that the cluster- 
ing strength, and hence the halo mass, is strongly correlated 
with the local density. The advantage of using SSFR instead 
of SFR is that SSFR measures the star formation efficiency 
directly, and hence it is easier to interpret the results. 

Interestingly, our results suggest that there are two popula- 
tions separated in the ro— SSFR diagram (Figure[5jl) by SSFR 
~ 2 x 10~ 9 yr -1 , which corresponds to the main sequence 
value reported at z ~ 2 (Dad di et ail 120071: iRodighier o et al.l 
120111) . For galaxies with SSFR < 2 x 10" 9 yr _1 (hereafter 
the low SSFR population), ro increases rapidly with decreas- 
ing SSFR, while ro increases mildly but significantly with in- 
creasing SSFR for galaxies with SSFR above the threshold 
(hereafter high SSFR population). 

4.2.1. Negative SSFR - ro relation: environment quenching of star 

formation 

The anticorrelation between ro and SSFR for the low SSFR 
population can be understood as an environmental effect simi- 
lar to that seen at lower redshif ts. For example, usi ng GALEX 
and SDSS samples at z < 0.3. iHeinis et all d2009) found that 
the galaxy clustering also declines strongly with SSFR. The 
explanation of those sBzK galaxies with very low SSFR in 
our sample are likely to be that they are falling into the denser 
environments where their star-formation activities are more 
effectively suppressed. Since their SSFR values are in be- 
tween the sBzK and pBzK galaxies, it suggests that these may 
be galaxies in transition between the main sequence and the 
quiescent population. And in fact, their clustering strength is 
indeed consistent with that of pBzK galaxies. 

On the other hand, since SSFR anticorrelates with stellar 
mass for sBzK as shown in the right panel of Figure|3] galaxies 
with lower SSFR tend to have higher M* and thus one might 
expect them to cluster more strongly, providing an alterna- 
tive explanation for the anticorrelation between ro and SSFR 

20 We note that although both lElbaz etal] UM%) and lCooper eTaTI <200l) 
claimed an inverted SFR-density relation at z ~ 1 , their results regarding the 
SSFR-density do not agree with each other. 



for low SSFR population. However, such an effect might not 
be dominant in the redshift range we are probing, as it ap- 
pears that the stellar mass dependence of SSFR in our sam- 
ple is weak and can not account for the strong correlation we 
see in the ro-SSFR relation. In order to further test whether 
the observed ro-SSFR relation is due to the combination of 
ro-stellar mass and SSFR-stellar mass correlations, we study 
the ro-SSFR relation in several stellar mass bins wherever we 
have enough statistics. As shown in Figure [7] we still see 
ro increases with decreasing SSFR for galaxies with SSFR 
below the critical value 2 x 10~ 9 yr" 1 , even in subsamples di- 
vided by stellar mass. We note that this result is robust against 
the incompleteness due to the missing of relatively quiescent 
galaxies (red z—K color) in the low stellar mass bins as men- 
tioned in §2 since they will only show up in the low SSFR end 
of the plot if they do exist. 

This implies that the ro-SSFR relation we see is not purely 
due to the stellar mass effect. Moreover, not only the sBzK 
galaxies with small SSFR, but the clustering strength of our 
pBzK galaxies is also found to be greater than that of the 
bulk of sBzK galaxies with similar stellar masses (see Fig- 
ure suggesting that quiescent galaxies preferentially re- 
side in denser environments. The negative ro-SSFR rela- 
tion we found for the low SSFR population is analogous to 
the monotonic decline in the mean overdensity with increas- 
in g SSFR that is seen at lower redshifts (e.g. see Figure 9 
of ICooper et al.ll2008l) . Furthermore, the higher clustering of 
pBzK galaxies relative to that of sBzK galaxies also mirrors 
the color-density relation found at lower redshifts in which 
the averaged over density is greater fo r galaxies with redder 
rest-frame colors (Cooper et afl 120061) . In other words, our 
results imply that some external processes that suppress the 
star formation activity in dense environments have started as 
early as z ~ 2. On the other hand, the positive correlation 
between ro and SFR as shown in Figure [5J: is likely due to 
the increasing population of massive star-forming galaxies in 
denser environment, simi lar to that found at z ~ 1 (Elbaz et al. 
120071: ICooperlFaIll2008l) . 

From Figure [5] it is also noted that the ro-SSFR has 
a steeper slope within the low SSFR range (SSFR < 2 x 
10~ 9 yr" 1 ), compared to the ro-M* relation, indicating that ro 
is more sensitive to SSFR than M*. This trend can be further 
illustrated in Figure|7] there exists a clear SSFR dependence 
of ro when splitting according to M*, while the variation of ro 
among different values is smaller for a fixed SSFR. Simi- 
lar effect was al so observed in a clu stering study of the SDSS 
sample done bv lHeinis et al.l (|2009), who claims that SSFR is 
a more sensitive probe of the halo mass than M*. Our result 
shows that such trend extends out to z ~ 2 as well. 

4.2.2. Positive SSFR — ro relation: environment triggering of star 

formation 

On the other hand, the trend we see for the high SSFR pop- 
ulation is not well-understood. According to their clustering 
strength, they have typical halo masses in the range 10 11 to 
10 13 M Q (see Tabled. However, Figure [3] shows that galaxies 
with the largest SSFRs tend to have relatively smaller stellar 
masses. Their stronger clustering therefore seems like a sur- 
prising deviation from the main ro - relation seen for the 
sBzK population as a whole (Figure^). Figure [7] shows that 
this positive correlation between ro and SSFR for high SSFR 
population is seen in all stellar mass bins. 

A possible explanation of our result is that the elevated 
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SSFR ("starbursts") is driven by mechanisms that are asso- 
ciated with more massive halos, even if the stellar masses 
of the starbursting galaxies are low. Galaxy interactions 
and mergers, for example, can trigger star formation activity 
(iBarton et al.1 ESobct lLambas et al.112003 1 iNikolic et all 120041: 
iLin et al.l 12007 : Ellison et al. 2008), and th ey are found to 
preferentially occur in denser envir onments dLin et alj|2010fc 
Ide Ravel et al]|2011t Uian et al.ll20T2h . In addition, the large- 
scale tidal field caused by grou ps/clusters can also induce star- 
bursts in interacting galaxies (Mart ig & Bournaudl2008l) . The 
working assumption is thus that the low stellar mass systems 
with enhanced SSFR are the ones located in denser environ- 
ments where interactions between galaxies are more common. 
If this is the case, one may question whether the high clus- 
tering amplitude can be attributed to the presence of close 
neighbors at small scales. However, we emphasize that the 
close pairs with angular separation less than 0.001 deg are ex- 
cluded when fitting the clustering amplitude, and therefore the 
high correlation length found for these high SSFR galaxies is 
directly linked to their large-scale environments. 

Infall shocking is another mechanism which is more ef- 
fective in massive halos. This may also compress the gas 
in galaxies, resulting in bursts of star formation. Such 
environment-driven enhanced star formation activity has 
also already been found in some c luster studies at z > 1 
dHilton et all 120101: iTran eTaT] 12010ft . We note that the en- 
hanced star formation activity in denser environments may 
not be easily seen in conventional environment studies which 
look for the quiescent fraction as a function of environment 
because the red (or blue) fraction may not change significantly 
even if the averaged star formation efficiency changes. 

Our finding concerning a dichotomy in the rg-SSFR rela- 
tion thus suggests that there are two opposite environmental 
effects influencing the star formation rate: quenching and trig- 
gering. While star formation can be reduced or quenched for 
galaxies located in denser environments, some other galax- 
ies residing in similar environments can have their star for- 
mation rates enhanced instead. Whether the mechanisms re- 
sponsible for these two effects are related or not is still not 
clear. One explanation is related to the well-accepted hypoth- 
esis that galaxy mergers, commonly found in denser environ- 
ments, can enhance star formation during the merger process, 
and then the remnants quickly become red and dead due to the 
lack of gas which is used up during the starbursting phase, or 
being blown out by the AGN activity. One way to test this pic- 
ture is to look for merger signatures of galaxies with elevated 
or suppressed SSFR. Recent study on the morphologies of 
starbursting galaxies at z ~ 2 already suggests that about 50% 
of these s ources are associated with interacting and merging 
galaxies (Karta ltepe et alj|201 1ft . It would be interesting to 
conduct a similar analysis for galaxies with suppressed SSFR 
as well. 

5. CONCLUSION 

We have taken advantage of deep HST/ACS data and 
CFHT/WIRCAM NIR data to identify z ~ 2 galaxies us- 
ing BzK color selection in the GOODS-N region. We have 
derived galaxy star-formation rates based on extinction cor- 
rected UV luminosity and calculated stellar masses with the 
Spitzer/IRAC 3.6pm and 4.5/^m photometry in order to study 
the clustering properties as a function of Zf-band magnitude, 
stellar mass (M*), star formation rate (SFR), and specific star 
formation rate (SSFR=SFR/M*). Our main conclusions are as 
follows: 



TABLE 1 

The values of ai and a 2 for 

THE EMPIRICAL FORMULA 
(EQ. [6]l TO COMPUTE THE 
STELLAR MASS. 



Redshift 


a\ 


a 2 


1.0<z< 1.5 


19.65 


-0.084 


1.5<z<2.0 


20.21 


0.069 


2.0<z<2.5 


20.79 


0.256 


2.5<z<3.0 


20.80 


0.890 



1. The UV-based SFR, corrected for extinction, scales with 
stellar mass as SFR oc M* Q , with a = 0.74 ±0.20 over the 
stellar mass range 10 9 Mq< M*< 10 11 M . The slope of 
the SSFR-M* is -0.26 ± 0.20, meaning a weak dependence 
of star formation rate efficiency on the stellar mass. 

2. The clustering amplitude of sBzK galaxies is a strong 
function of K-b&nd magnitude and stellar mass, increasing for 
more massive and brighter galaxies. 

3. Highly star-forming galaxies are more strongly clustered 
than galaxies with low star formation rates, which is in line 
with the correlation between SFR and stellar mass. 

4. For the first time, we find that the correlation length 
reaches a minimum at SSFR of 2 x 10~ 9 yr _1 , the typical value 
for the "main sequence" of star-forming galaxies at z ~ 2. The 
correlation length is larger for galaxies with both smaller and 
larger SSFRs. Such a dichotomy holds even at fixed stellar 
mass. Our results suggest that environment has two effects: 
quenching and inducing the star formation activities. Stronger 
clustering for galaxies with relatively low SSFR implies that 
environment has started playing a role in quenching star for- 
mation at z ~ 2, while another environment effect, galaxy 
interactions and mergers, might explain the elevated SSFRs 
("starbursts") in more massive halos (denser environment). 

5. Passive galaxies {pBzKs) are more strongly clustered 
than sBzK galaxies at a given stellar mass, suggesting that 
the color-density relation is in place at z ~ 2. The correlation 
length ro of pBzK galaxies is measured to be 24.6±7.0 h~ l 
Mpc, which is larger than that of sBzK galaxies but is similar 
to that of sBzK galaxies with the lowest SSFRs. 

6. Our results suggest that current Durham semi-analytical 
models of galaxy formation appear to underestimate the SFR 
and predict larger scatter in SFR for star-forming galaxies of 
a given stellar mass. As a result, the predicted ro-SFR and 
ro-SSFR relations cannot be directly compared to the ob- 
served results, although the K and dependence of clus- 
tering is in better agreement with observations. 

One caveat in our analysis, however, lies in the small sam- 
ple size that prevents us from a more comprehensive HOD 
analysis. The ongoing Spitzer Extended Deep Survey (SEDS; 
PI: G. Fazio), combined with existing multi-wavelength data 
in several extra-galactic fields, will provide larger samples 
with better determined photometric redshifts and stellar mass 
measurements. This survey will thus allow improved mod- 
eling of the relationships between halo mass and the galaxy 
properties through the combined clustering and abundance 
analysis. 
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TABLE 2 

Clustering properties and corresponding dark matter halo masses of sBzK and pBzK galaxies in GOODS-N. 



Sample Cut 


z 


number 


Aj a) 






1 OA 
DM 




M m 
DM 








2 21 


ftfi 

OU 


n ni 497n -i- n nnzts^zi 

U.UlfZ l\j HI U.UU4Jjt 


17 4 zh 4 9 


5 2 ± 8 


/7 o _i_ 9 q\ v 

^ / .Z ZTZ Z.7^ A 




(5.0 ±2.3) 


x 


10 13 


(4. 9 ± 2.3" 


x 10 13 


91 S <^ 99 n 


9 1 % 

Z.l J 


1 83 


n nnszi^i 4-n nnmszi 

U.UUJtLll IL U.UUlUo't 


in 3 4- 9 s 
iu.j zrz z.j 


q a _i_ n a 
j. j zrz u. j 


^ i . y zrz u. \j) a 


in 13 


(1.1 ±0.4) 


x 


10 13 


(1.0 ±0.4) 


x 10 13 


11 (\ ^ K «f 99 S 


9 1 S 

L. 1 J 


9QS 

Z 7 J 


n fYnnnR 4- n nnn^ss 


7 zl 4- 1 Q 

/ .*+■ ZTZ 1 .y 


9 4 4- n ^ 
z.h zrz u. j 


(-7 i _|_ 7 c\ v 
i^/.i zrzz.o^ a 


in 12 

1U 


(3.4± 1.5) 


x 


10 12 


(3.0 ± 1.5) 


x 10 12 


99 S ^- / 9"! fj 

ZZ.J <-.. A v, Z.J.\J 


2 24 


441 


n nn9n9s 4- n nnnzioi 

U.UUZUZJ JT_ u.uuutyi 


S R 4- 1 S 
j.o zrz i .j 


2 ± 2 


/a 7 _i_ i ft\ v 
i^j.z zrz i a 


in 12 


(1.3 ± 0.8) 


x 


10 12 

1U 


(1.1 ±0.7) 


x 10 12 


9^ A «< V «*-* 1% S 


9 1Q 

z. i y 


709 
/ uz 


n nn9 s i zi 4- n nnno s/^ 

U.UUZJ 1*+ JZ U.UuUZjO 


f\ (\ 4- 1 S 
O.O ZC 1 . J 


9 9 4- n 1 

Z.Z ZTZ U. 1 


rs n 4- 1 c\\ v 

ip.U ZTZ 1 .UJ A 


in 12 

1U 


(2.2 ±0.5) 


x 


10 12 
1U 


(1.9±0.5) 


x lO 12 


9-2 s ^- <^ 9zl 


2 12 


i n9Q 

1UZ7 


n nm 1 s i 4- n nnm ^s 

U . UU 1 1 J 1 jZ <j . UUU 1 J 


at.-},-] n 

t. j zrz i .u 


1 s 4- n 1 
i.j zrz u. i 


1^7. J ZTZ D. L r) A 


in 11 

1U 


(3.4 ± 1.4) 


x 


IQll 


(2.1 ± 1.1) 


x 10" 




1 78 
l . / o 


1688 


n nn 1 1 S7 4- n nnn i nzi 


zl S 4- 1 n 
t. j zrz i .u 


1 S 4- n 1 
i . j zrz u. i 


i^i . j zrz u. j ) a 


in 12 

1U 


(4.7 ± 1.2) 


x 


in 11 

1U 


(2.9 ±0.9) 


x 10" 


q s <- 1no(M f\A^ \ <r inn 


2 13 


1201 


n nni 817 4- n nnm ^s 

U.UU 1 O 1 / JZ U.vUv 1 J J 


5.6 zh 1.2 


1 9 ± 1 


(2.7 ± 0.4) x 


10 12 


(1.1 ±0.2) 


x 


10 12 


(8.9± 1.7) 


x 10" 


lOO/ InaiM f\A^. \ 1fl ^ 


9 I'X 


/ jZ 


n nn99Qi 4-n nnn9ss 
u.uuzzyi zc u.uuuzjj 


ft -3 _(_ 1 /I 

o. j zrz i .t 


9 i 4- n 1 

Z.l zTZ yJ. 1 


1 zrz u.yj A 


in 12 


(1.8 ± 0.5) 


x 


10 12 


(1.5 ±0.4) 


x 10 12 


10S ^ JnofM f\/f^\ / 1 I fl 

iu.j <^ /(^^yw^/iviQj <, ii. u 


9 9£ 




n nnssss 4- n nnnft^fi 

U.UUjooo ZC U.UUUOjO 


1 n £ 4- 9 zl 

lu.o zrz z.'f 


j.j zrz u.z 


(1 \ 4- n A\ v 

l^Z. 1 ZC U.^f _J A 


in 13 


(1.2±0.2) 


x 


10 13 


(1.1 ±0.2) 


x 10 13 


Jn-ofM f\A^ \ ^> 1 1 n 


2 27 


1 07 


n nno^^9 4- n nn9Q6i 

u.uuvooz zn u.uuzyo 1 


ij.7 zrz j.y 


/I q _|_ n 7 
t. j zrz u. / 


(■zL ^ 4- 1 8^ v 
l_t. j zrz i.oj a 


in 13 


(2.7 ± 1.3) 


X 


10 13 


(2.6 ± 1.3) 


x 10 13 


1 n <" SFR/f'M^vr - ^'! / s n 
i.u v. oriv^iviQ yi ^ j.u 


1 84 


OJ 1 


n nm S97 4- n nnni 

U.UU1JZ / X U.UUUlOU 


S Q 4- 1 9 
j . j zrz i .z 


i 7 _l n 1 
i . / zrz u. i 


i^z.j zrz u.o^ a 


in 12 

1 U 


(9.1 ±2.7) 


X 


10" 


(6.8 ±2.5) 


x 10" 


5 < SFR/fM^vr~h < 10 


2.06 


1013 


0.001392 ± 0.000162 


4.9zh 1.1 


1.7±0.1 


(1.5 ±0.4) x 


10 12 


(5.9 ±1.8) 


X 


10" 


(4.3 ±1.6) 


x 10" 


lO<SFR/(M yr i )<3O 


2.13 


1243 


0.001814 ±0.000133 


5.6 ±1.2 


1.9±0.1 


(2.7 ±0.4) x 


10 12 


(1.1 ±0.2) 


X 


10 12 


(8.9±1.8) 


x 10" 


3()<SFR/(M yr- i )<6O 


2.27 


484 


0.002571 ± 0.000437 


6.6 ±1.6 


2.2 ±0.2 


(5.0 ± 1.6) x 


10 12 


(2.3 ±0.8) 


X 


10 12 


(2.0 ±0.8) 


x 10 12 


6O<SFR/(M yr" 1 )< 100 


2.33 


203 


0.005680 ± 0.001059 


10.2 ±2.5 


3.3 ±0.3 


(1.9 ±0.6) x 


10 13 


(1.1 ±0.4) 


X 


10 13 


(1.0 ±0.4) 


x 10 13 


SFRAMoyr 1 ) > 100.0 


2.21 


218 


0.009798 zh 0.002099 


14.1 ±3.5 


4.4 ±0.5 


(4.5 ± 1.3) x 


10 13 


(2.8 ±0.9) 


X 


10 13 


(2.8 ±0.9) 


x 10 13 


2 x 10" 10 < SSFR/yr" 1 < 4.5 x 10" 10 


2.33 


105 


0.018360 ± 0.002629 


19.6 ±4.6 


5.8±0.4 


(9.6 ± 1.6) x 


10 13 


(6.8 ±1.4) 


X 


10 13 


(6.8 ±1.4) 


x 10 13 


4.5 x 10" 10 < SSFR/yr" 1 < 7.5 x 10" 10 


2.22 


146 


0.005615 dz 0.001473 


10.3 ±2.7 


3.3 ±0.4 


(1.9 ±0.8) x 


10 13 


(1.1 ±0.5) 


X 


10 13 


(1.0 ±0.5) 


x 10 13 


7.5 x 10" 10 < SSFR/yr" 1 < 1.5 x 10" 9 


2.13 


708 


0.001666 zh 0.000267 


5.3 ±1.3 


1.8±0.1 


(2.2 ±0.8) x 


10 12 


(8.9 ±3.6) 


X 


10" 


(6.9 ±3.2) 


x 10" 


1.5 x 10" 9 < SSFR/yr" 1 < 3.5 x 10" 9 


2.13 


1883 


0.001141 zh 0.000090 


4.3 ±1.0 


1.5 ±0.1 


(9.1 ± 1.8) x 


10" 


(3.3 ±0.7) 


X 


10" 


(2.1 ±0.6) 


x 10" 


3.5 x 10" 9 < SSFR/yr" 1 < 7 x 10" 9 


1.79 


847 


0.003525 zh 0.000308 


8.4±1.9 


2.6±0.1 


(1.3 ±0.2) x 


10 13 


(6.2 ±1.1) 


X 


10 12 


(5.6 ±1.0) 


x 10 12 


SSFR/yr" 1 > 7 x 10" 9 


2.38 


301 


0.005726 zh 0.000639 


10.2 ±2.3 


3.3 ±0.2 


(1.9±0.3)x 


10 13 


(1.0 ±0.2) 


X 


10 13 


(1.0 ±0.2) 


x 1() 13 


pBzK's 


1.99 


44 


0.025320 zh 0.008292 


24.6 ±7.0 


7.1 ±1.2 


(1.5 ±0.6) x 10 14 


(1.2 ±0.5) 


x 10 14 


(1.2±0.5) 


x 10 14 



NOTE. — All values listed above are for the sBzK galaxies except for the last row, which is for the pBzK galaxies. (fl> Correlation amplitude at 1 degree; ib) Correlation length in h 1 Mpc; (c l Bias ; ld) 
Halo mass in Mq inferred directly from the Mo & White 1 2002) formalism for the value of ro ; <e> Minimum halo mass in Mq assuming ro is the effective clustering strength averaged over halos above 
a certain threshold; ^ Minimum halo mass in Mq assuming that ro is the effective clustering strength averaged over halos above a certain threshold and weighted by the applied HOD, as described in 
the text. 
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